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Abstract. We study the long-time behavior of the 3-dimensional repulsive nonlin- 
ear Hartree equation with an external attractive Coulomb potential —Z/\x\, which 
is a nonlinear model for the quantum dynamics of an atom. We show that, after a 
sufficiently long time, the average number of electrons in any finite ball is always 
smaller than AZ (respectively 1Z in the radial case). This is a time-dependent gen- 
eralization of a celebrated result by E. H. Lieb on the maximum negative ionization 
of atoms in the stationary case. Our proof involves a novel positive commutator 
argument (based on the cubic weig ht |x| 3 ) and our findings are reminiscent of the 
RAGE theorem. 

In addition, we prove a similar universal bound on the local kinetic energy. In 
particular, our main result means that, in a weak sense, any solution is attracted to 
a bounded set in the energy space, whatever the size of the initial datum. Moreover, 
we extend our main result to Hartree-Fock theory and to the linear many-body 
Schrodinger equation for atoms. 

© 2012 by the authors. This paper may be reproduced, in its entirety, for non- 
commercial purposes. 



1. Introduction and Main Result 

The rigorous attempt to answering the question "How many electrons can a nucleus 
bind?" has received substantial attention in the literature over the last decades [4T1I421 
EMEU 54, 57, 58, 62, 63 . So lar, the question was only addressed in a time-independent 
setting, that is, the absence of bound states was shown when the number of electrons 
in the atom is too large. In the present paper we shall rigorously formulate and provide 
an answer to a similar question in the time-dependent setting: " How many electrons 
can a nucleus keep in its neighborhood for a long time?". 

Our main purpose is therefore the rigorous understanding of the long-time behavior 
of atoms. We shall prove for instance that, in the Hartree approximation, a nucleus of 
charge Z cannot bind in a time-averaged sense more than 4Z electrons (respectively 2Z 
in the radial case). In particular, we will recover some of the known time- independent 
results (non-existence of bound states) by different arguments, though. One key in- 
gredient in our paper turns out to be a new commutator estimate leading to a novel 
monotonicity formula, which may be of independent interest for both linear and non- 
linear Schrodinger equations. 

As a model for the quantum dynamics of an atom, let us first consider the time- 
dependent nonlinear Hartree equation with an external Coulomb potential: 

i x) = ( -A - + M 2 * 7^7 ] u(t, x), 



(1.1) 



dt v ' ' \ M \x\ 
u(0,x) = u (x) £ H 1 ^ 3 ). 

Here u(t, x) describes the quantum state of the electrons (which are treated as bosons 
for simplicity) in an atom [3011311110] • The terms in the parenthesis are respectively the 
kinetic energy operator of the electrons, the electrostatic attractive interaction with 
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the nucleus of charge Z, and the mutual repulsion between the electrons themselves (in 
units such that m = 2 and h = e = 1). The total number of electrons in the system is 
a conserved quantity, which is given by 

/ \u(t,x)\ 2 dx = \u a (x)\ 2 dx =: N. 

JM 3 J R 3 

In physical applications, the number N is an integer but it is convenient to allow 
any positive real number here. Note that, in Section [4] below, we will also consider the 
physically more accurate Hartree-Fock model as well as the full many-body Schrodinger 
equation describing atoms. But for the time being, we deal with the Hartree equation. 

The nonlinear equation (|l.ip and many variations thereof have been studied exten- 
sively in the literature. The existence of a unique strong global-in-time solution to 
(jl.ip with an initial datum uq £ H 1 (M 3 ) goes back to Chadam and Glassey [13]. Their 
argument is based on a fixed point argument combined with the conservation of the 
Hartree energy defined by 

E z (u):=[ Wui^dx-Z [ l -^fdx + U [ \<f<f d X dy. (1.2) 
J»3 J R3 \x\ 2 J R3 J R 3 \x - y\ 

In fact, the global well-posedness result for (II. ip can be extended to initial data in 
L 2 (R 3 ); see for instance [111I5S] . However, in what follows, we will always assume that 
Mo lies in the energy space iJ^R 3 ) so that its corresponding energy is well-defined. 

When Z ^ 0, the solution u(t) to the Hartree equation (11 . If) exhibits a purely 
dispersive behavior, which has been studied by many authors. Here, some works were 
devoted to the understanding of the dispersive effects for any initial datum 17,;22 l2"51l3"2l 
1331155] , whereas several others dealt with the construction of (modified) scattering fM\ - 

In this paper, we are interested in the physically more relevant case when Z > 
holds, which corresponds to having an external attractive long-range potential due to 
the presence of a positively charged atomic nucleus. The electrons can (and will) now 
be bound by the nucleus, and the problem to understand the long-time behavior of 
solutions is much more delicate. For instance, it was already noticed by Chadam and 
Glassey in [13l Thm. 4.1] that the solution u(t) cannot tend to zero in L°°(R 3 ) as 
t — >• oo for negative energies &z{uo) < 0, which can occur if Z > holds. 

When Z > 0, there exists nonlinear bound states, which are solutions of Equa- 
tion (jl.lj) taking the simple form u(x)e~ ltx , where u £ i/ 1 (M 3 ) solves the nonlinear 
eigenvalue equation 

For any fixed < N ^ Z, it is known that the equation (|1.3[) has infinitely many solu- 
tions such that J K 3 \u\ = N. Moreover, there is a unique positive solution, which mini- 
mizes the Hartree energy (|1.2p [3,43. subject to N fixed, and the other (sign-changing) 
solutions can be constructed by min-max methods 44 , 64 , 68 . The interpretation of 
the condition < N ^ Z is that the atom is neutral (if N = Z) or positively ionized (if 
N < Z). In this situation, it is not energetically favorable to send a positive fraction of 
L 2 -mass /i > 0, say, to spatial infinity, since the remaining charge is Z — (N — fi) > 
positive and thus attractive far away from the origin. A more precise mathematical 
statement is that the Palais-Smale sequences with a bounded Morse index cannot ex- 
hibit a lack of compactness when N ^ Z, and this implies the existence of infinitely 
many critical points [51 1231145] . 

It is known that there are bound states in the case of negative ionization, i. e. when 
N > Z holds. By [201 Thm. 7.19] (see also 0[S]), there is a minimizer of the Hartree 
functional for N slightly larger than Z. However, it is physically clear that there should 
not be any bound state when TV is too large compared to Z , because a given nucleus is 
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not expected to bind too many electrons compared to its nuclear charge. In [Tl HOETj . 
it was proved that there exists a universal critical constant 1 < 7 C < 2 such that the 
equation (|1.3p has no solution for N > ^y c Z but has at least one for N ^ J C Z- That j c 
is independent of Z follows from a simple scaling argument. 

Let us now collect some basic facts about the set of solutions of the time-independent 
problem (|1.3p . For any u E H 1 (M 3 ), the self-adjoint operator 



Z ■ ,2 # J_ 

x\ ' \x\ 



-A- A 1 '-' a 



has no positive eigenvalue, by the Kato-Agmon-Simon theorem [501 Thm. XIII. 58]. 
This shows that necessarily A in Equation (11.31) . Furthermore, we can derive an 
upper bound on ||Vu||l2 which only depends on Z as follows. If u € i? 1 (R 3 ) solves 
(|1.3[) . then by taking then the scalar product with u we find that 

f \Vu(x)\ 2 dx^Z f ^ff^ dx^Z \\Vu\\ L2 \\u\\ L2 . 

Here we have used the inequality 

L dx * So (I I ^ 2 + h 1 1 Vm|2 ) = M ^ l|Vu W 3 > 

which follows from the value of the hydrogen ground state energy, inf Spec (—A/2 — 
= — z 2 /2. We conclude that any solution u e H 1 {R 3 ) to (|1.3|) must satisfy the 

bound 



|Vu| 2 < lc Z\ 

Recalling that L 3 \u\ 2 ^ C Z holds, we conclude that the set of all stationary states 

s/ z ■■= [u E i^QR 3 ) : u solves fL3) for some A sC oj (1.4) 

is bounded in 7J 1 (IR 3 ). Elementary arguments show that s^z is weakly compact in 
i/ 1 (IR 3 ). But we note that the set srfz is not compact in the strong i/ 1 -topology. 

Supported by physical reasoning and rigorous results in linear scattering theory 
about asymptotic completeness (see Remark [1] below) , it is common belief for infinite- 
dimensional Hamiltonian systems such as (jl.ip that any of its solutions should behave 
for large times as a superposition of one or several states getting closer to the global 
attractor s$Zi plus a dispersive part. This is what has already been shown for Z ^ 0, 
in which case srfz = {0}. Not much is known in this direction for nonlinear Schrddinger 
equations [65, 66 , and solving this problem (a. k. a. soliton resolution) constitutes a 
major mathematical challenge. For the Hartree equation (jl.ip studied in this paper, 
the situation is even less clear because of possible modified scattering due to the long- 
range effects of the Coulomb potential. We can, however, formulate a simpler (but 
weaker) conjecture as follows. 

Conjecture 1 (The global attractor). Let u(t) be the unique solution to the Hartree 
equation (|1.1[) for some uq £ H 1 (M 3 ). Take any sequence of times t n — > oo such that 
u(t n ) — * it* weakly in i? 1 (IR 3 ). Then u* € s^z- 

Remark 1 (The many-body Schrddinger case). Let us recall that the Hartree equa- 
tion (jl.ip is a nonlinear approximation of the linear many-body Schrddinger equation 
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On the contrary to the Hartree case where we can allow N — J„ 3 \u\ 2 to take any 
positive real value, the number N of electrons must of course be an integer for (|1.5|) . 
The Hartree equation (jTTTJ) is obtained by constraining the solution ^S(t) to stay on the 
manifold of product states of the form ^(t, X\, ...,xn) — tp(t,xi) x ••• x ip(t,XN) and 
using the Dirac-Frenkel principle. Then u(t) — \^Nip(t) solves (11.11) . Let us remark 
that (|1.5[) can be rewritten after a simple rescaling as follows 




Thus the limit of large N —¥ oo with N/Z fixed corresponds to the usual mean-field 
limit. In this regime, Hartree's theory is known to properly describe (bosonic) atoms, 
both for ground states J6}j and in the time- dependent case J4> 19L See also J21, 561 for a 
review on mean-field limits and the Hartree approximation. 

The many-body equation (|1.5[) looks complicated, but it has the advantage of being 
linear. In particular, the RAGE theorem tells us that the only possible non-zero weak 
limits ofi&(t) when t — > oo are bounds states of the Hamiltonian H(N) in the parenthe- 
sis fI \\18[{ST\ \ 5i $ . This is not a very precise description of the solution for large times 
because if some particles stay close to the nucleus while other escape to infinity, we will 
always get ^f(t) — 1 weakly in H 1 (M. 3N ), see '391. However, asymptotic completeness 
is known to hold for the linear evolution equation (|1.5[) . This exactly says that any so- 
lution ty(t) is, in an appropriate sense, a superposition of bound states of the operators 
H(k) with I ^ k ^ N and of scattering states l761|ff7[[5P) /. Because of the behavior of 
the underlying many-body system, it is reasonable to believe that the same should be 
true for the Hartree equation p. II) . 

A somewhat weaker property that would follow from Conjectured] (at least for (II. 7|l ) 
is that for large times, the local mass of any solution has to be smaller than j c Z. 

Conjecture 2 (Asymptotic number of electrons and kinetic energy). Let u(t) be the 
unique solution to the Hartree equation (|1.1[) for some Uq G //^(R 3 ). Then 



and 




for all r > 0. 

The upper bound jcZ 3 is certainly not optimal here. In physical terms, the con- 
jecture says that whatever the number of electrons we start with (and whatever their 
kinetic energy), we will always end up with at most j c Z electrons having a universally 
bounded total kinetic energy. The other electrons have to scatter because the attraction 
of the nucleus with positive charge Z is not strong enough to keep all the electrons in 
its neighborhood. It could be that proving the weaker Conjecture [2] is not much easier 
than proving the stronger Conjecture [TJ We actually have very little information on j c . 

In this paper, we are interested in Conjecture [5J We will prove a time-averaged 
version of (jTTTJ) , with j c replaced by 2 in the radial case, and by 4 in the general case. 
Our main result is as follows. 

Theorem 1 (Long-time behavior of atoms in Hartree theory). Suppose Z > 0, let uq be 
an arbitrary initial datum in H 1 ^ 3 ), and denote by u(t) the unique solution of (|1. II) . 
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Then, for any R > 0, we have the following estimate 

1 f T , f , Ht,x)\ 2 Arr 3 2y r KNR 2 

— dt dx 1 v , J) „ ^ 4Z + — + — (1.9 

T J J R 3 1 + \x\ 2 /R 2 R ZT 

N:= f Kl 2 



with 



and 



(1.10) 



K := sup f \Vu(t)\ 2 ^Z 2 N + 2 \Vu \\l m3) + N 3 \\Vu \ 
In particular, we have 

limsupi f dt I dx \u(t,x)\ 2 < 4Z (1.11) 

T^OC 2 Jq ,/|x|<r 

for every r > 0. Similarly, we have the following estimate on the local kinetic energy 

i -i T dti dx) vu ^:<(^ + 2 ^ + '4)^r T dt r dx ^ 2 



Tj J m (1 + N/i?) 2 ^ V 4 i? ff/Tij 7 K 3 l + \x\ 2 /R 2 

2R*/Ks/N 



T (1-12) 



and therefore 



limsupi / dt /" |Vw(t,x)| 2 s$ Z 3 

/or every r > 0. 

J/ £/ie initial datum uq = uq(\x\) is radial, then u(t) is radial for all times and the 
same estimate (ll.9[) holds true with 4Z replaced by 2Z . Similarly, the estimate (|1.13p 
holds true with Z 3 replaced by Z 3 /2. 

Note that we do not exactly get that the limiting mass is ^ 4Z for large times, but 
we only know it in the sense of time averages of the form (f)r = f dt. Such a 

statement is reminiscent of the celebrated RAGE theorem [TJ[TH1[5T1[52] for linear time 
evolutions generated by self-adjoint operators. The constants in the error terms of (II. 9[) 
and (|1.12p are probably not optimal at all, but they are displayed here to emphasize 
that our method can provide simple and explicit bounds. However, we have not tried 
to optimize these constants too much. 

In the radial case, we are able to get the same numerical value of 2 as the best known 
estimate on j c . However, we use a virial-type argument that seems to be quite different 
from Lieb's celebrated proof in [4T] in the stationary case (which, for radial solutions, 
goes back to Benguria [7]). In particular, our approach provides an alternative proof 
of the fact that j c < 2 in the stationary radial case. 

Strategy of the Proof. Now, we explain the main ideas used in the proof of Theo- 
rcm[TJ To this end, we start by quickly recalling Lieb's proof [41] that j c < 2 holds. His 
idea is to take the scalar product of the stationary Hartree equation (|1.3p with |a;|u(a;), 
leading to the estimate 

u M-a) + (-a)\z\\_ zn+ [ r M±MM^l dxdy<0 , 

2 / its i 8 3 2|x-y| 

using that A ^ holds. To conclude, it suffices to notice that we have 



(1.13) 



|x|(-A) + (-A)H = 1/2 



w) 1x11/2 >0 > 
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by Hardy's inequality, and that 

M + M ^ x 
\x-y\ 

by the triangle inequality. Combining these estimates, we obtain that — ZN+N 2 /2 < 0, 
which implies the bound N < 2Z for the stationary problem (jl.3p . (Note that the 
inequality is strict, since there is no optimizer in Hardy's inequality.) 

In view of Lieb's argument for the stationary problem (|1.3[) . it appears to be a viable 
strategy in the time-dependent setting to consider the quantity M(t) — J \x\ \u(t, x)\ 2 dx 
(or some spatially localized version thereof). Indeed, if we take the second time deriv- 
ative of M(t), we are (formally) led to the well-known Morawetz-Lin-Strauss estimate 
for nonlinear Schrodinger equations (NLS), which has been proved of enormous value in 
the setting of NLS with purely repulsive interactions. However, due to presence of the 
attractive term —Z/\x\ with Z > in the Hartree equation the use of the classical 

Morawetz-Lin-Strauss bounds does not yield any dispersive information about u(t, x), 
even in the case when TV is large compared to Z. 

In our situation, it turns out that it is more natural to study the time evolution of 
the third moment M(t) — J \x\ 3 \u(t, x)\ 2 dx. If we compute its second time derivative, 
we obtain 

I ah L |a;|3|u(t ' X)|2 dX = Jt {u{t)lAu{t)) 

= (u(t),i[-A,A]u(t)) + (u(t),i[V u ,A]u(t)) 

with A := -i(V ■ x \x\ + \x\x ■ V) and V u = -Z\x\^ + \u\ 2 * This is the same 

as multiplying the time-dependent equation (jl.l[) by Au(t) and taking the imaginary 
part. Our key observation is the positivity of the commutator 

i[-A,A] = -~[A,[A,\x\ 3 }] >0 (1-14) 
(see also (|2.6I) below), combined with the fact that 

(u(t),i[V u ,A}u(t)) =-2 1 \x\x- \JV u {x)\u{t,x)\ 2 

(\x\x-\y\y) ■ \u(t, x)\ 2 \u(t, y)\ 2 dx dy - 2ZN 

> kN 2 ~ 2ZN, 

where k = 1 if u{t) is radial and k = 1/2 otherwise (see Lemma [3] below) . Hence, when 
N > 2Z/k, we deduce the lower bound 

lis f M 3 \u(t,x)\ 2 dx=^-{u{t),Au(t))^N(KN-2Z)>0. (1.15) 
3 at z J R 3 dt 

Therefore the quantity J R3 |a;| 3 |u(t,a;)| 2 dx grows at least like t 2 for large t and in 
particular (u(t),Au(t)} is a monotone increasing quantity. This growth is a strong 
indication that some dispersion takes place and some particles have to escape to infinity. 
(A regularized version of the previous estimate will indeed show this claim for any Ab- 
solution.) Note also that, in the time-independent case when u is a nonlinear bound 
state (and hence the left side in (|1.15l) must be zero), this is also a new proof of 
Lieb's inequality j c < 2 in the radial setting, since k = 1 holds under this symmetry 
assumption. 

Let us generally remark that virial or positive commutator arguments are very com- 
mon in the literature [U1GI5]. When |a;| 3 is replaced by \x\ this leads to the famous 
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Morawetz inequalities [27J as already mentioned, whereas the case of \x\ 2 gives the 
virial identity used by Glassey in |29] to prove finite-time blowup for NLS. In a recent 
work [66], Tao advocated the use of \x\ A for some nonlinear Schrodinger equations in 
dimension d 7, in order to get a universal bound on the mass of the solution. We are 
not aware of any use of the multiplier \x\ 3 in the literature. 

In fact, using the cubic weight |cc| 3 is rather natural from a dimensional point of view 
in our situation: If the potential term [V^,j4] should be O(l), then the virial function 
must behave like the third power of a length to compensate the Laplacian and the 
Coulomb potential. 

For the proof of our main result, we will in fact derive a whole class of double 
commutator estimates of the same kind as (|1.14p . which we think is of independent 
interest too. In particular, we will show in (|2.6[) below that, in any dimension d 1, 
we have the commutator bound 

- [A, [A, |*|"]] + d - 4)(d - /3)|*r 4 , (1.16) 

provided that (3 ^ max(l,4 — d). Note that the right side is ^ when (3 ^ d. In 
spite of the fact that (|1.16[) turns out to be equivalent to a general version of Hardy's 
inequality, we have not found it explicitly written (let alone systematically treated) in 
the literature. Notice that the bound (|1.16D contains the usual inequalities for f3 — 1, 2, 
as well as Tao's estimate for (3 = 4. In the present application, we shall use (|1.16|) 
in dimension d — 3 with j3 = 3, or rather a regularized version thereof. However, let 
us remark that the positivity of this commutator does not directly follow as in the 
"classical" cases when fj = 1, 2. To wit this, we note that, for d = f3 = 3, a calculation 
(which will be detailed below) yields the identity 

- [A, [A, |x| 3 ]] = -AA|x| 3 - V • (Hess N3 )V 
24 

= -p- 12V- [\x\(l + w x u;J)]V, 

where lu x = denotes the unit vector in direction x € M 3 . Obviously, the first term 
on the right side is negative definite. Nevertheless, when combined with the second 
term, the generalized Hardy's inequality (see (|2.7p below) shows that we indeed have 
that the whole right-hand side is non negative, and hence the estimate (|1.16[) follows 
in the particular case d = f3 = 3. 

Ultimately, we are interested in general 7J 1 -solutions u(t) without imposing any 
spatial weight condition. Therefore, the strategy of proving Theorem [1] explained above 
needs to be further refined. In particular, the desired bound p.9|) on a ball of radius 
R cannot be obtained by only looking at the second derivative of the third moment as 
we have just explained. Our method to extend (|1.9|) to any H 1 -valued solution u(t) is 
to replace the function |x| 3 by a radial function /i?(|a;|) which behaves like |a;| 3 on the 
ball of radius R and like \x\ at infinity. This will imply that Af R = — i[A, /r] defines a 
bounded operator on L 2 (R 3 ). Furthermore, we will need to derive a sufficiently good 
lower bound on the double commutator —[A, [A, /#]] in order to imitate the previous 
argument on the ball only. In Section[2] we explain how to do this for a general function 
/. Finally, the bound (|1.12[) on the local kinetic energy is itself obtained by considering 
another virial function gn which behaves like \x\ 2 on the ball of radius R and like \x\ 
at infinity. The complete proof of Theorem [T] is given in Section [3] 

Extensions: Hartree-Fock and Many-Body Schrodinger Theory. In physi- 
cal reality, electrons are fermions, which means that the many-body wave function 
= ty(t, x\, xn) in (|1.5[) must be antisymmetric with respect to exchanges of its 
spatial variables x\, ...,xn- The Hartree state ip(t,xi) ■ ■ ■ t/j(t, xn) is symmetric and it 
is therefore not allowed for physical electrons. This is why one speaks about bosonic 
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atoms. The simplest product-like antisymmetric wave function is a Hartree-Fock state 
sometimes also called a Slater determinant 

i f!(t,x 1 ,...,xjf) = -y= ^2 e(a)u 1 {t,x a{ x))---u N {t,x^ N) ), (uj,u k ) L2 = S jk . 

In Section 14.11 below, we extend Theorem [T] to the corresponding time-dependent 
Hartree-Fock equations; see Theorem [5] for a precise statement. Finally, we also con- 
sider the full many-body Schrodinger equation p.5[) in Section 14.21 below, where our 
findings are summarized in Theorem [3] 



2. Estimating the Commutator — [A, [A, /(x)]] 
Throughout this section, we use the convenient notation 

p := -iV, 

and in particular we have p 2 = —A in what follows. In this section, we investigate how 
to get lower bounds for a double commutator of the form — [p 2 , [p 2 ,f(x)]\ in general 
space dimensions d ^ 1. Such a double commutator always arises when computing the 
second derivative of the expectation value of f{x), in a non-relativistic system based 
on the Laplacian. We always assume that / is smooth enough (possibly only outside of 
the origin), such that the double commutator can be at least properly interpreted as a 
quadratic form on C™(R d ) or on C^(R d \ {0}). 

Our starting point is the well-known formula for the double commutator, which 
follows from a tedious but simple calculation: 

- [p 2 , [p 2 J(x)}] =-(AAf)(x)+ip-(Eessf(x))p. (2.1) 

Since the Hessian of / appears on the right side, it is natural to restrict to convex 
functions /. Then the second term is non-negative in the sense of operators. One can 
use this term to control the bi-Laplacian of / by resorting to Hardy's trick, which is 
based on writing 

p ■ (Hess f(x)) p= (p + iF(x)) ■ (Hess f(xj) (p - iF(x)) 

+ i(p- (Hess f(x))F(x) - F(x) • (Hess f{x))p) 

- F{x) ■ (Ressf(x))F(x) 

^ div (Kessf(x)F(x)) - F(x) ■ (Hess f(x))F(x) (2.2) 

for any sufficiently smooth real vector field F : M 3 — > R 3 . Here we have only used 
that {jp + iF(x)) ■ (Hess /(a;)) (p — iF(x)) ^ holds, which simply follows from the 
assumed convexity Hess/(x) ^ and the self-adjointness {p + iF(x))* = p — iF(x). 
For dimensional reasons, it is natural to take F of the form F(x) = ax\x\~ 2 with some 
constant a € M. We thus obtain the lower bound 

-[p 2 , lp 2 J(x)}] =i(p + iaj^j ■ (Hess /(a)) (p - ia-^j 

+ 4a div (^Hess/(x)^) - ^a gXgggMf _ {AAf){x) 
> 4a div (Vss/t*)^) - 4q 2 * T ( H ^/(*))* _ (AA/)(x) (2 3) 

for a sufficiently smooth convex function / and any a € E. Note that by using Hardy's 
trick we are able to obtain a lower bound which does not contain the differential operator 
p. Our estimate only involves a multiplication operator. By varying a, we can try to 
make the negative part of this function as small as possible. 
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Let us now restrict ourselves to a radial function /(|a;|) and use the notation r = \x\ 
and uj x := x/\x\ for simplicity. Some tedious calculations show that 



Hess f(\x\) = (1 - ^ w J)ffl + ^/"(r), 



r 



div ((He SS / W )i)^v(m.,)^ + (,- 2) ™ 
x T (Ressf(x))x _ f"(r) 



r-2 



\x\* r 

Moreover, we recall the formula for the Bi-Laplacian of a radial function: 

AA/(|ar|) = f^(r) + 2(d - + (d-l)(d~ 3)^^ - (d - l)(d - 3)^. 

Therefore we can rewrite the equality in (12.31) for a radial function / as follows: 

-/W( r ) + 4fa- d -^ /(8) W 



2 (d-l)(d-3)\ /"(r) 



2 

4 I a(d - 2) - a 2 - 

+ (d-l)(d-3)^. (2.4) 

The operator on the first line is ^ when x i— > /(|x|) is convex. In dimension <i = 3, 
we already get a simple estimate. 

Lemma 1 (A lower bound for d = 3). Let f : [0, oo) — > K 6e a convex non- decreasing 
function such that x i— ► /^ 4 )(|a;|) € L 1 1 oc (R 3 ). TTien we Aave 

-[/, [ P 2 ,/(N)]] =4( P + J ^) • jV^J)^ + ^J/>)) 

-/ (4) (N) 

>-/ (4) (kl) (2-5) 
in the sense of quadratic forms on C^°(M. 3 ). 

Proof. Take a = 1 in J23J). □ 

Coming back to (|2.4|) and taking now the convex function /(|x|) = \x\° with (3 1, 
we obtain the following general result. 

Lemma 2 (Estimate on — [p 2 , [p 2 , \x\P]]). For all (3 max(l,4 — d), we have 

- [p 2 , [p 2 , \xf}} 2 f3((3 + d - 4)(d - p) \xf~\ (2.6) 

in the sense of quadratic forms on C^°(R d ) (resp. on C™(M. d \ {0}) if 13 = 4- d). The 
right side of (|2.6p is non negative for max(l, 4 — d) j3 ^ d. 

Proof. Take f{r) = r* 3 in (|2 .4[) and optimize with respect to a (the optimum is a = 
(j3 + d — 4)/2). We need /3 ^ 1 to make sure that / is non-decreasing and convex, and 
(3 > 4 — d to ensure that all the terms are in L 1 1 oc (E rf ). For (3 = 4 — d ^ 1, the right side 
of (|2.4[) vanishes and the bound stays correct by a simple limit argument. We remark 
that, in the borderline case (3 = 4 — d, there is a positive 5-measure occurring at the 
origin x = 0, which we do not see when using functions of C^°(M. d \ {0}). □ 
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Remark 2. Note also the special formula — [p 2 , [p 2 , \x\ 2 ]} — 8p 2 ^ valid in any 
dimension d 1, which immediately follows from (|2.1j) . For d ^ 3 and (3 = 2, the 
lower bound given in Lemma&is then a direct consequence of Hardy's inequality 4p 2 ^ 
(d — 2) 2 \x\~ 2 . In fact, we shall see below that the bound in Lemma\^is equivalent to a 
generalized version of Hardy 's inequality. 

We conclude this section with some general observations as follows. First, we note 
that Lemma [2] gives a non negative lower bound in (|2.6[) in dimension d = 2 for the 
choice (3 = 2 only. In higher dimensions d ^ 3, the right side is non negative for any 
1 ^ (3 ^ d. When (3 = 4, we get the simple lower bound 

-[p 2 ,[p 2 ,\x\ 4 ]} ^ 4d(d-4), ford>4, 

which was used for the first time by Tao in [66] . 

As we have seen, the bound (|2.6[) is equivalent to the operator inequality 







(p + 


-t) - ( 




= r 13 . 


This can ; 



^ 



as 



\xf- 2 [\PtVu(x)\ 2 + (/3-l) 



r, / \ (3 + d-A . , 

uj x ■ Vu(x) H — j — u(x) 

2\x\ 



dx > 0, 



where P^- = 1 — w^wj is the projection on the two-dimensional space orthogonal to ui x . 
Saying that the second term is non negative is equivalent, for (3 > 1, to the (generalized) 
Hardy inequality 



(2.7) 



f \xf- 2 \u x -Vu(x)\ 2 dx^ + d 4)2 / |x|^ 4 Kx)| 2 dx. 

JR d 4 J Rd 

Hence we see that (|2.6p is nothing else but a reformulation of Hardy's inequality (|2.7 
Remark 3 (Fractional Laplacians). Using the integral representation 



g sin(7r#) 



■ s^ds, for < 9 < 1, 



we can easily transpose most of our estimates to fractional powers \p\ 20 = (— A) e and 
(p) 26 = (|f>| 2 + l) 9 with 9 G (0,1). For instance, for the pseudo-relativistic kinetic 
energy operator \J p 2 + 1, we have, at least formally, 



^i+p 2 ,yi+p 2 ,f{x)\ 



; ds 



y/idt> 



OTWTWTWTT) (- [p2 > [p2j{x)]] ) (i 



i 



+ P 2 + s){i+p 2 + ty 



In particular we find 



VT+pi, [Vi + ^,\xf] 



> o 



for max(l,4 — d) ^ (3 ^ d. For a general convex radial function f and in d = 3 
dimensions, we obtain the estimate 



v/TT^, [VT+^,/(|x| 



i 



(4) 



with /l denoting the positive part of f^ . 
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3. Proof of Theorem [T] 

In this section, we provide the proof of our main result given by Theorem [TJ We 
always assume that the initial datum uq is smooth and decays fast enough, such that 
our calculations are justified. As we will see below, our estimates only involve the 
H 1 (M. 3 ) norm of uo, and thus the general case can therefore be obtained by a simple 
limiting argument, which we do not detail here. 

Step 1. The Virial Identity. Consider a smooth radial convex function /. We define 
the corresponding virial operator 



A f := p ■ V/ + V/ -p = p ■ oj x f(\x\) + /'(MK • p. 
Using Formula (I2.5[) of the previous section, we get 



(3.1) 



dt 



[u(t),A f u(t)) =4 



P^Vu(t,x)\ 2 dx 



+ 4/ /"(M) 



uj x ■ Vit(t, a;) + 



u(t, x) 



dx 



with 



f {i X\x\)\u{t,x)\ 2 dx-2 / f(\x\)\u(t,x)\ 2 u x -VV u (t ) x)dx, 

(3.2) 



V u {t,x) = — n + \u{t)\ 2 *\x\- 1 :=— n + W u {t,x). 
\x\ \x\ 



The first potential term is just 

/'(M)|u(i, x)\ 2 oj x ■ V (y-^j dx = ~ 2Z 
The second potential term can be expressed as follows 
-2 f f{\x\)\u{t,x)\ 2 u x -VW v {t,x)dx 



CM. 

S.3 \A 2 



u{t,x)f 



= 2 



x-y 



/'(NIK , 

3 \x-y\ A 
(f(|x|K-f(M)o; y )-(x-y) 
\x - y\ 3 

where in the last line we have just exchanged the role of x and y. Inserting in (|3.2p . we 
arrive at the following expression 



u(t,x)\ 2 \u(t,y)\ 2 dxdy 

u(t, x)\ 2 \u(t, y)\ 2 dx dy 



dt 



[u(t),A f u(t)) =4 



K3 \x\ 

n\x\) 



P^Vu(t,x)\ 2 dx 



U! x ■ Vu(t, x) 



u(t,x) 


2 

dx 


\x\ 


2Z [ , 


(N) 



(f'(\x\) Ux -f'(\y\) Uv )-(x-y) 
\x - y\ 3 



'dx 



u(t, x)\ 2 \u(t,y)\ 2 dxdy. 



(3.3) 



For dimensional reasons, it is natural to take /(|a;|) = |x| 3 /3. The following lemma 
allows to deal with the last potential term in this special case. 
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Lemma 3 (Lower bound on the nonlinear term for /(r) = r 3 /3). We have 

(\x\ 2 uj x - \y\ 2 uj y ) -{x-y) 1 

\x-y\* "2 { -> 

for all x ^ y € M 3 . In the radial case we have 

-f (|l|V ', ly|2 f (X ' y) ^ du y = 1 (3.5) 

where f s2 duj x = (Att)^ 1 j s2 doj x denotes the (normalized) angular integration. 

Proof. We compute 

{\x\ 2 uj x - \y\ 2 UJ y ) - jx-y) _ r 3 + s 3 - (r 2 s + s 2 r)uj x ■ uj y 
\x — y\ 3 [r 2 + s 2 — 2rsoj x ■ u y ) 3 / 2 

_1 + u 3 - (u + u 2 )9 
~ (l + u 2 - 2u6) 3 / 2 

with x = ruj x , y = su) y , u :— min(r, s) max(r, s)~ l G [0, 1] and 6 :— lo x ■ u) y G [—1, 1]. 
Differentiating with respect to 8, we find 

d (1 + u 3 - (u + u 2 )6\ u(l + u){u 2 + (1 - u)(2 - it)) - 6»u 2 (l + u) 



d6> V (1 + « 2 ~ 2u6») 3 / 2 J (1 + u 2 - 2u6fl 2 

We have u 2 + (1 — u)(2 — it) = u + 2(u — l) 2 tt and therefore the numerator is non 
negative for u > and 9 6 [—1,1]. We conclude that the minimum is attained for 
8 = —1. The value is 

1 + u 3 + u + u 2 _ 2u 1 

(1 + u) 3 =1 ~ (1 + u) 2 ^ 2 

where the minimum is attained for u = 1. All in all, we find that 

(\x\ 2 w x - \y\ 2 uj y ) -(x-y) I 
\x-y\ 3 ' V 

as was claimed. In the radial case we find by explicit integration 

f f (\x\ 2 cj x \y\ 2 uj y ) .( X -y) J_ f 1 + u 3 (u + u 2 )0 M = 

Js*Js* \^-y? V 2J-1 (1 + u 2 - 2u6) 3 / 2 

This concludes the proof of Lemma [3] □ 



For /(r) = r 3 /3, the previous estimates gives 



^-(u(t), A f u{t)) = 4 f \x\ \P±Vu(t, x)\ 2 dx + 8 f \x\ 



U! x ■ Vu(t, x) + 



u(t, x) 



dx 



+ kN 2 - 2ZN (3.6) 



where k = 1 in the radial case and k = 1/2 otherwise. If u is a stationary state, then 
the left side is independent of t and this is a new proof that N < 4Z (resp. 2Z) for 
bound states. Equation (|3.6p is a new monotonicity formula for the Coulombic Hartree 
equation, when N 4Z (resp. iV 2Z in the radial case). 
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Step 2. The Localized Virial Estimate. We now use a localized virial estimate, 
which means that we choose a virial function which behaves like |x| 3 /3 on a ball of 
radius R and like |x| at infinity. We will take fa of the form 

f R {\x\)=R 3 f(\x\/R) 

for 

f(r) = r — arctanr, (3-7) 

which we have chosen to have 

f\r) = -^ = l-—^. (3.8) 
w 1 + r 2 1 + r 2 y ' 

Clearly, the first derivative /' is non-decreasing and positive. Hence x H> /(|x|) is a 
convex function on IR 3 . The following lemma gathers some important properties of /, 
which are the 'localized' equivalent of Lemma [3] above. 

Lemma 4 (The virial function /). Let f be as in (|3.7[) . We have 

(f'(\x\)u> x -f'(\y\)u y )-(x-y) l f(\ X \) f(\y\) 

\x-y\* " 2 |x| 2 |y| 2 

for all x ^ y e M 3 . In the radial case, we get 

(f'(\x\)u> x -f'(\y\)u y ).(x-y) 

\x^yf X " 

/'(max(|x|,M)) _ 1 f'(\x\) f'(\y\) 

max(|x|,|y|) 2 1 + max(|x| 2 , \y\ 2 ) ' |x| 2 \y\ 2 ' K ' ' 

Proof. Like in Lemma |3l we write 

(f'(\x\)u x - f'(\y\)w v ) • (x - y) _ rf'(r) + sf(s) - 9(sf(r) + rf'(s)) 



(3.11) 



|x - y\ 3 (r 2 + s 2 - 2rs9) 3 / 2 

with r = |x|, s = \y\ and 9 = u> x ■ u> y £ [— 1, 1]. Differentiating with respect to 6, we find 

r(2r 2 - s 2 )/'(r) + s(2s 2 - r 2 )f'(s) - 9rs(sf'{r) + rf'js)) 
(r 2 + s 2 - 2rs6») 5 /2 ' 

Since /' > 0, the numerator is positive for 8 ^ 9 C and negative for 9 ^ 9 C . Regardless 
whether 9 C E [—1, 1] or not, the minimum of the function in (|3.11[) is attained at 9 = ±1. 
For 9 = — 1, we find 

/'(r) + /'(a) _ r 2 + .s 2 + 2r 2 s 2 
(r + s) 2 _ (1 + r 2 )(l + s 2 )(r + s) 2 ' 

Now we remark that 

r 2 + s 2 + 2r 2 s 2 _ 1 ( (r - s) 2 + 4r 2 s 2 1 



(r + s) 2 2 2(r + s) 2 " 2 : 

and therefore 



(r + s) 2 " 2(l + r 2 )(l + s 2 )' 
For = 1, we find 

\f'(r)-f'(s)\ = |r 2 (l + .s 2 )-s 2 (l + r 2 )| = r + s f'{r)f{s) 
(r — s) 2 (1 + r 2 )(l + s 2 )(r — s) 2 \r — s\ r 2 s 2 

We have, with u = min(r, s) max(r, s) -1 , 

r + s 1 + u 2u 

= 1 + > 1. 



r — s 1 — u 1 — u 
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We conclude that 

(/' (MR - f'(\y\)w y ) ■ (x - y) > 1 f'(\x\) f(\y\) 
\x-y\* ' 2 \x\ 2 \y\ 2 ' 

for all x 7^ y € K 3 , as was stated. 

In the radial case we have to compute the integral over the angle explicitly. We use 
the notation r< := min(r, s) and r> :— max(r, s), and we get 

(f(\x\)u x -f'(\y\)u v )-(x-y) 



is 2 J s 2 



\x-y\* - x ^y 



1 f 1 rf(r) + sf(s)-9(sf'(r)+rf'(s))_ de 



_! (r 2 + s 2 - 2rs6») 3 / 2 

r/'(r)+*/'(«)l T d9 1 



2 7_ x (l + u 2 -2ufl) 3 / 2 
(r/'00 + */'(r)) 1 f 1 



d0- 



r?, 2}_ x (1 + u 2 - 2m6») 3 / 2 

r/'frO + s/'OO 1 (rf( S ) + sf(r)) u 

r 3 1 — u 2 r3 1 — M 2 

1 



r>(r 2 - r 2 ) 
1 



(r>/'(r>) + r</'(r<) - (r</'(r>) + r > /'(r < ))r < /r > ) 



, ^r(r>/'(r>)-^/'(r>)/r>) 



r>(r 2 - r 2 ) 
_ f(r>) 

Note also that the previous calculation is valid for an arbitrary radial differentiable 
function /, and not just the specific / chosen above. The proof of Lemma U is now 
complete. □ 

We apply p.3[) for f R — R 3 f(-/R) with / given by (I3.7[) . We get the expression 
j t (u(t),A fR u(t)) 



m f RI ' { \ X \ lR) \PtVu{t,x)\ 2 dx + AR [ f\\x\/R) 
Jr3 F Jr 3 



, . u(t, x) 2 
oj x ■ \/u(t,x) H j — j — ace 



-4/ / (4) (k|/ J R)Kt^)| 2 dx-2z/ fl2/ ^J /fl) | M (*^)l 2 ^ 

•« J R 3 \X\ 2 

+ i? / / ^ : - '- \u{t,x)\ \u(t,y)\ dxdy. 

Jm Jm F - 2/1 

(3.12) 

We now define the localized mass by 



M R (t):= if!-»\u(t,x)\ J dx= ——^\u(t,x)\"dx. (3.13) 
J R 3 |x| 2 J R3 1 + i? 2 |a;| 2 

Using (|3.9p (resp. (|3.10[) in the radial case), we get the lower bound 

d 1 f 

-{u(t),A fR u{t)) > -- / f^(\x\/R)\u(t,x)\ 2 dx- 2ZM R (t) + nM R {tf (3.14) 
with k = 1 in the radial case and K = 1/2 otherwise. Finally we remark that 

~/ (4, M = ^T^7 > - 24 (T^ 1 < r > '» "IT? - ^ 
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(the best numerical constant is 1.33 instead of 3) and we get our final lower bound 

j t (u{t),A fR u(t)) > - ( 2Z + |) M R (t)+KM R (t) 2 . (3.15) 

To conclude our proof of (|1.9[) . we average (13.151) over a time interval [0,T] and use 
Jensen's inequality 

1^ M R (t) 2 dt> (1^ M R (t)dt\ , 

to get 

(u(T),A fR u(T))~(u(0),A fn u(0)) 



T 

> K [h lv Mli{t)d *) ~ {2Z + 3/R) ('r lo 



M R {t)dt\. (3.16) 



Note that 

Kt*(t),^/«u(t))| = |(u(t),(p-v/H(N) + v/fl(N)-p)u(t) 

sC 2y/Ky/N\f R \ Laa = 2VkVnR 2 

since sup r>0 f'{f) = 1, and where we recall that K — sup t ||Vu(t)| i 2. In summary, we 
conclude that 

k (jf J o M R (t) dt\ - (2Z + 3/R) (± jT M R (t) dt \ < 4 Vi^V ^. 
Using y/1 + u ^ 1 + u/2, this implies 



1 f T 

- J M R (t)dt^ 



2Z + 3/R 2Z + 3/R / 16kVKNR 2 



2k 2k V (2Z + 3/i?) 2 T 



2Z 3 AVKNR 2 2Z 3 2VKNR 2 
^~ + R + (2Z + 3/R)T ^~ + R + ZT ' 

which ends the proof of (|1.9[) . 

Remark 4. 0«r proof works exactly the same for a more general time average based 
on a positive function fi such that J Q /i — 1 and fx' is a bounded Borel measure. More 
precisely, we have the estimate 



00 



,(t/T) dt f H t, X )\ 2 ^2_Z + 3 + VKNtf 



T 7 R 3 1 + \x\ 2 /R 2 k R ZT 

For instance, one could take fi{t) = e _t . 

Step 3. Estimate on the Local Kinetic Energy. We show here that the kinetic 
energy also has a universal upper bound in average, on any ball of radius R. This time, 
we use a localized virial identity based on the function 

g R (\x\)=R 2 g(\x\/R) 

which behaves like \x\ 2 on B R and like \x\ at infinity. More precisely, we take 

g(r) = r - log(l + r) (3.17) 

which is such that 

r 1 
9'(r) = — - = 1 



1 + r 1 + r 

Clearly g' is positive and non-decreasing, therefore x 1— > ff(|x|) is convex on R 3 . 
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We use the lower bound (|2.4p with a = and we get, by the same calculations as 
before, 

^(u(t),A gR u(t)) =*f ( ^^ \P x L Vu(t,x)\ 2 +g''(\x\/R)\u }x -Vu(t,x)\^dx 



2^/ MI^ W M)| 2 ^ (3 . 18) 



We denote by 



\Vu{t,x)\ 



2 



K R (t):= g"(\x\/R)\Vu(t,x)\ 2 dx = y ' n dx 

Jr3 Jr3 (1 + i? ^xl) 2 

the local kinetic energy. Since x i— >• g(|ic|) is convex, then (Vg(x) — V g(y)) ■ [x — y) ^ 
for all x, y £ M 3 . Also, we notice that 

" U (1 + r) 2 ^ 1 + r r 

Finally, we compute 
and 



WM = _5_<J_-./'( r ) 

(1 + r) 4 ^ 1 + r 2 



ff W (0 = - " r = 

So we arrive at the estimate 



j t (u(t),A gB u(t)) > AK R (t) - jgM R (t) - 2zJ R9 ' [ ^J R) \u{t,x)\ 2 dx. (3.19) 
In order to control the negative term, we use again the trick of Hardy: 



0^ / g"{\x\/R) \Vu(t,x) + auj x u(t,x)\ 2 dx 
in 3 



a 2 



g"{\x\/R)\Vu{t,x)\ 2 dx- a / div(u x g" (\x\/R)) \u(t, x)\ 2 dx 

g"(\x\/R)\u(t,x)\ 2 dx 
g"(\x\/R)\X7u{t,x)\ 2 dx + a 2 I g" (\x\/R)\u(t, x)\ z dx 



2a 



*' {l *\ /R) \u{t,x)\ 2 dx-?- [ g^(\x\/R)\u(t,x)\ 2 dx. 
\x\ R J M 3 



Therefore, using that ~5 (3) (r) = 2(1 + r)~ 3 2(1 + r 2 )" 1 = 2/'(r)r" 2 and that 
g"{r) = (1 + r)~ 2 s$ (1 + T 2 )- 1 = f'{r)r- 2 , we find 

9"i\A/R)s„.u _m 2j _^ 1 ^ m , ^ , 1 



\u{t,x)\ z dx^—K R (t)+\- + -)M R {t) 
R 3 |x| 2a \ 2 R 



DYNAMICAL IONIZATION BOUNDS FOR ATOMS 



17 



Coming back to the negative term in (|3.18[) , we write 



R9 'W/ R \u(t, X )\> dx 



fi3 \X\ 2 

1 



\x\(l + \x\/R) 



\u(t, x)\ 2 dx 



I 1 f 1 

\u(t, x)\ 2 dx + — / j — \u{t, x)\ 2 dx 



\x\{l + \x\/R) 2 ' v ' n Rj m (l + \x\/R) 



2a " w \2 ii 
Inserting in (|3.19[) gives 

j t {u{t),A gR u{t)) > (^-^K R {t)-Z^oc+^ + ^)M l M). (8.20) 

Taking a = Z/2 leads to 

A<u(t),A te tt(t))>2iir Ji (t)-zf| + | + ^) (8.21i 

To conclude our proof, we average over t in an interval [0, T] using that 
\{u{t),A gR u{t))\ < 2>/KVF|0k| io . = 2RVKVN 

and we get 

which concludes the proof of (|1.12j) . 

Step 4. Estimate on K. We end the proof of Theorem[T]by estimating the maximal 
value K of the kinetic energy of u(t), in terms of |uo|a-i, using the conservation of 
energy. 



Lemma 5 (Kinetic energy estimate). We have 



ll Vu (<)li,2(R3) < Z 2 \\uq\\ 2 L 2 {r3) + ||Vuo||^ 2(R3) + - \\u f L2(R3) ||Vuo|| i2(HS) (3.23) 
for all t G M. 

Proof of Lemma\^ By conservation of energy and mass, we find 

A Z\ \ 1. 



£z(u ) = £. z (u) ^ (u, ( -— - — ) u) + t;I Vu Il 2 (r 3 ) 



>-4" / \u \ 2 + h\Vuf 



2 

since -A/2 - Z^^ 1 ^ -Z 2 /2 (hydrogen atom). Next, for x G R 3 and u e H 1 
we note the bound 



|z - 2/1 
which gives us 



rfy sC mm ( ~ f \u\ 2 + f \Vu\ 2 ) = |u| L2(R3) ||Vu|| i2 



(3.24) 



2 (R3) ' 2 H"0||L2(R3) || V «U|| L 2( R 3) • 

Hence, 

|V«li2 (H 3) «S ^ I«0li2(R8) + 2 |VUo|l2 (R3) + |W0li2 (R 3) || Vli || L 2 (R 3 ) . 

□ 
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This concludes the proof of Theorem [TJ □ 



4. Extensions: Hartree-Fock and Many-Body Schrodinger Theories 

4.1. Hartree-Fock Theory. The Hartree-Fock equations describe the nonlinear evo- 
lution of a wave function taking the form of a Slater determinant, i. e., 

*(*) = ~y= ^2 sgn(cr) -ui(t,x CT(1 )) •••ttw(i,a; .(jv)) 

where the functions ui, un model the states of the N electrons. The physical fact 
that electrons are fermions is expressed in the Pauli principle given by the orthonor- 
mality condition 

(uj,Uk) L 2 = Sjk- 

The Hartree-Fock equations 10 12 Hi] form a system of N coupled nonlinear equations 
similar to (jl.ip 

>,—Uj = H uUj , 

N N (4.1) 

H u v = -Av - Z\x\~ x v + \ u k\ 2 * \x\~ lv ~ * \ x \~ lu k- 

fe=i fe=i 

One simple way to write the same equation is to introduce the one-body density matrix 

N 

7(*) : = \ u k)( u k\ 

k=l 

which is the orthogonal projection onto the space spanned by the functions ui, ...,un- 
Then (|4.1I) is equivalent to the so-called von Neumann equation 

d_ 

'at' 



*^77 = [H~f,j], 

r (4-2) 
H^v = -Av - Z\x\ 1 v + p l{t) *\x\ 1 v- I j(t,x,y)v(y)dy. 



Here Pj{x) := "/(x,x) is the density associated with the matrix 7. The time-dependent 
equation (|4.2[) makes in fact sense for any trace-class operator 7 such that 

sC 7 s$ 1 and Tr( 7 ) = N, 

which corresponds to generalized Hartree-Fock states [2]. Note that the infinite rank 
case rank 7 = +00 is also allowed here. We refer to [TUJ[T2] for the proof of global 
well-posedness for (|4.2[) with initial data such that Tr(l — A)7 < +00. 

The following result is the equivalent of Theorem [T] in the Hartree-Fock case. 

Theorem 2 (Long-time behavior of atoms in Hartree-Fock theory). Suppose Z > 
and let 70 be an arbitrary initial datum such that 

Tr(l - A) 7o < 00. 

Denote by j(t) the unique solution of (14. 2[) . Then we have the following estimate 

1 f T , f , PtmO) 3 2VKNR 2 , . 

' dt dx 71 <k 4Z + 1 + - + — (4.3) 



T J 7 R 3 1 + \x\ 2 /R 2 R ZT 

with 

N := Tr( 7o ) 

and 



K := supTr(-A)7(i) < Z 2 N + 2Tr(-A) 7o + N 3 y/Tr(-A)-y Q . (4.4) 
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In particular, we have 

rT 



limsupi f dt f dx p 1 {t){x) < AZ + 1 (4.5) 

T->oo T J J\ x \<r 



' \x | 

/or every r > 0. Similarly, we have the following estimate on the local kinetic energy 
lf T dt[ dx- ^ 



ry v (i+H/i?) 2 



where t 7 (x) = — Ylk=i {®kldk){x, x) is the density of kinetic energy, and therefore 

1 f T f Z 2 
limsup— / dt dxT l{t) (x) ^ — (4Z + 1) (4.7) 

/or every r > 0. 

// i/ie initial datum 70 is radial in the sense that 

10 {@x,&y)= 7o (x,y), Vx,yeR 3 , V^e50(3), 

£/ien 7(i) is radial for all times and the same estimates (|4.3[) and (|4.7D /10/d irwe wit/i 
4Z + 1 replaced by 2Z + 1 . 

The proof of Theorem [5] is very similar to that of Theorem [TJ the main new difficulty 
is the control of the exchange term. Thus we only explain how to deal with it. 

Sketch of the Proof of Theorem^ First, we consider a sufficiently smooth radial func- 
tion / = /(|ir|). (Below we will take / = //?,, the same as in the proof of Theorem [T]) 
Differentiating with respect to t, we find 

j t Tr(A n )=iTr([H 1 ,A f h) 
„2 rji 



Tr([p 2 ,[p 2 ,/]] 7 ) +zTr([T/ 7 ,A / ] 7 ) - iTr([X^A f ] 7 ) (4. 



where V 1 — ~Z\x\ 1 + |x| 1 * /? 7 and X 1 is the exchange term defined by 



(X y u)(x)= [ f^Au^dy. 
Jm? \x-y\ 

Note that 

i[Vy,A f ] = -2V/- W 7 

is a function (that is, a multiplication operator). Analogous to the Hartree case, we 
thus obtain 



iTr([V^,A f ] 1 )=-2 / Pl (x)Vf{x)-VV 1 {x)dx 
= — 2Z ( ^-—Tn— p^it, x) dx 



(V/Qe) - V/(y)) ■ (s - y) 

j \i Pi 0*0 p 1 [y) dx dy. 

\x - y\ a 

The exchange term is controlled using the following fact. 

Lemma 6 (Exchange term). Let Tr(l — A) 7 < +00 and suppose f : M. d —> K satisfies 
V/ G L°°(R d ). Then we have 

iTr([X 7 ,A fh )=[ f (V/W-V/(y))-(x-.) 2dx% 
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Proof. The proof is an explicit computation: 

iTr([X^,A f ] 7 ) =iTr([X 7 , (p- (V/) + (V/)-p)] 7 ) 

X 7 (x, y)V y ■ (V/)(y)7(y, x) da; dy 

l{y, x)V x ■ (V/) [x)X 1 (x, y) dx dy 

x 1 (x ) y)(Vf)(y) ■ v y j(y,x)dxdy 

l{y,x){V f)(x) ■V x X 1 {x,y)dxdy. 
Integrating by parts for the first two terms we find 

iTr([X 7 ,A f }~/)=- f f _j(y,x)(Vf)(y)-VyX^x,y)dxdy 

X 7 {x,y)(Vf)(x) -Vx7(y,x)dxdy 

x j( x ,y)(^f)(y) ■^ y iiy,x)dxdy 

l(y,x)(Vf){x) ■ \J x X 1 {x,y)dxdy. 
Now we use that 

^yXjix, y) = ^ V,7(a, y) + l(x, y)V„ * ■ 

and we exchange x and y in the second and fourth integral. The final result is 
iTr([X^A f } 7 ) 

\j(x,y)\ 2 ( (V/)(y) • V ^/ ^^ + (V/)(ar) • V.j-^) dxdy 

V F-J/I \x-y\J 

(V/(x) - V/(y)) ■ (x - y) 

i 13 7(^2/) a^ay- 

\x - y\ 6 

□ 

Inserting this in (|4.8p gives the following value for the derivative of the expectation 
value of Af. 

^-Tr(A fl ) = -Tr(b 2 ,[p 2 ,/]] 7 ) - 2Z f ^M^i,*) dx 

(V/(x)- V/( y ))-(x-y) (p7(a;)p7(2/) _ |7(X)J/)|2) ^ (4 1Q) 
p 2/1 

Since / is convex, we have the operator bound 

-b 2 ,[ P 2 ,/]]W (4) (M), 

which gives 

-Tr([p 2 ,[p 2 ,/]] 7 ) > -Tr(/( 4 ) 7 ) = - / (\x\) p l(t) {x) dx 

because of 7 ^ 0. Thus we can argue exactly as in the Hartree case. We start by taking 
fa given by p.7[) and define the local mass by 

M R {t):= [ tML Pl{t) { x )dx. 

JR 3 \ x \ 
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Then we use the bound (|3.9[) . that is 



(v/flW - v/«(y)) ■ - v) i gVWg) g/Wg) 

|aj-y|3 ^2 \x\ 2 |y| 2 ' 

as well as the fact that p 1 {x)p 1 {y) ^ \~i(x,y)\ 2 for a. e. x,y € K 3 (by the Cauchy- 
Schwarz inequality and the eigenfunction expansion for 7.) This gives 

(Vfe " ) -, VM ; )H '-» l w ,kw - hferio ** 

with /ifl := R 2 f'(\x\/R)\x\~ 2 . Since < 7 < 1 and < h R < 1, we have h R ^h R < 
(7ir) 2 ^ ^ii and therefore 

Tr(/i R 7/ifl7) < Tr(^ 7 ) = M R {t). 

We conclude that 

2 

The additional term is responsible for the change of 4Z into 4Z + 1. In the radial case, 
we use (|3.10[) instead and we get rid of the factor of 1/2 on the left-hand side. The rest 
of the proof is exactly the same as in the Hartree case. □ 

4.2. Many-Body Schrodinger Equation. Our method also applies to the linear 
many-body Schrodinger equation 

i^-*(t) = H(N,Z)*(t), 

U(0) = * GH 1 (( K3 ) Ar )' 
of which the Hartree and Hartree- Fock models are nonlinear approximations. 

The Hamiltonian H(N, Z) is self-adjoint and bounded from below on L 2 ((M. 3 ) N ) with 
domain H 2 ((R 3 ) N ) and quadratic form domain iJ 1 ((R 3 ) JV ). Of particular interest are 
its restrictions to the symmetric (a. k. a. bosonic) and antisymmetric (a. k. a. fermionic) 
subspaces. These are also self-adjoint operators, denoted respectively by H S (N, Z) and 
H a (N, Z). In cither of these two subspaces, the essential spectrum of H a / S (N, Z) is a 
half line [E„/ a (JV, Z), 00) where 

Z a/s (N, Z) = inf Spcc(H a/s (N - 1, Z)), 

by the HVZ-Theorem [15150] . It is known that there are no positive eigenvalues [2U] , but 
there might be embedded eigenvalues in [E / s (iV, Z),0]. There exists a critical number 
of particles N°, S (Z) such that H a j s [N,Z) has no eigenvalues below E / s (iV, Z) for 
N > N^ /s {Z), see [511571155] . For bosons, it is known that 

hm S±-L = 7c 

where j c ~ 1.21 ^ j c is the largest number of electrons that ground states can have in 
Hartree theory [5j|6l[6T|. For fermions, it was proved in [42] that 

lim *P=1. 

Z-¥<X> Z 
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The best bound valid for all N goes back to Lieb [JT] and it holds both for bosons 
and fermions: N°,(Z) < 2Z + 1. For fermions, it was recently improved to N£(Z) < 

1.22 Z + 3Z 1 / 3 by Nam [49]. 

All the previous authors seem to have only studied when the Hamiltonian H a i s (N, Z) 
ceases to have eigenvalues below its essential spectrum. The question of the existence 
of embedded eigenvalues in [E a / S (N, Z),0] does not seem to have been addressed so 
far. But this is a relevant problem in the context of the time-dependent equation. Our 
method allows us to prove that there are no eigenvalue at all when N ^ 4Z + 1. 

Theorem 3 (Linear many-body Schrodinger equation). The Hamiltonian H(N, Z) has 
no eigenvalue when N AZ + 1 . 

Here we do not distinguish between the different particle statistics. Thus our result 
applies to all of L 2 ((R 3 ) N ) and it deals with all possible symmetries. We, however, 
conjecture that the largest N such that H a / S (N,Z) can have eigenvalues behaves like 
N c a/s {Z) for large Z. 

Proof. Let * g H 2 ((R 3 ) N ) be an eigenfunction of H (N, Z) and let f R ( | x | ) = R 3 f ( | x \ /R) 
be as in p.7[) . Then we write 

0=(*,i (H(N,Z)J2(A fR ) Xj -J2( A fnUH(N,Z)J *j 



N N I / 

= E <*. W (^.M*) - 2 E (*. *M*j) • v., - R + - 2 E ^-^| 

j=l j=l \ \ fe/j 1 J 1 



^ / y- (V/fl(^j) ~ V/fl(x fc )) ■ (xj - x fc ) . ^ 



Using p.9p . we get 




i?^(N)#%(W) 



|x fe | 2 



HI 



"] I 2 



p^(x)dx\ —— j — ,2 py(x)dx. 



In the last line we have used Jensen's inequality as well as the fact that f'(r)/r 2 — 
1/(1 + r 2 ) < 1. Passing to the limit as R — !• oo gives N <AZ + 1. □ 

Since H(N,Z) has no eigenvalue when N ^ AZ + 1, it follows from the known 
existence of scattering and the asymptotic completeness 16,37,59] that any solution 
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^f(t) of the time-dependent equation (|4.11[) behaves (in an appropriate sense) as a 
superposition of bound states of H (k, Z) with k < 4Z + 1 plus a scattering part. In 
particular, it is possible to prove that 

limsup / py(t, x) dx ^ 4Z + 1. 

By using argument in the proof of Theorem [3] and following step by step the method 
of Section [3J one can get a simple proof of the weaker result 

limsup — f dt J p^{t, x) dx 4Z + 1. 

T^oo T Jq J\ x \^r 
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